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Abstract 

We fit the reduced cross section for deep-inelastic electron scattering data to a three parameter 
In 2 s fit, A+(5\d 2 {s/ sq), where s = x)+m 2 , and Q 2 is the virtuality of the exchanged photon. 

Over a wide range in Q 2 ( 0.11 < Q 2 < 1200 GeV 2 ) all of the fits satisfy the logarithmic energy 
dependence of the Froissart bound. We can use these results to extrapolate to very large energies 
and hence to very small values of Bjorken x — well beyond the range accessible experimentally. 
As Q 2 — > oo, the structure function F^x, Q 2 ) exhibits Bjorken scaling, within experimental errors. 
We obtain new constraints on the behavior of quark and antiquark distribution functions at small x. 

PACS numbers: 13.60.Hb, 12.38.-t, 12.38.Qk 
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Introduction. Inclusive deep-inelastic lepton scattering (DIS) has played a seminal role 
in particle and nuclear physics, notably for its early manifestation of Bjorken scaling 
and, soon thereafter, the logarithmic scaling violations that are a hallmark of perturbative 
quantum chromodynamics. From the structure functions of DIS, one also obtains crucial 
constraints on the parton distribution functions essential for predictions of short- distance 
hard-scattering phenomena at very high energies, whether at hadron collider facilities or 
in ultra high energy cosmic ray interactions. The DIS lepton-nucleon cross sections may 
also be analyzed in complementary fashion as hadronic Vp scattering cross sections, where 
V = 7*, a virtual photon in the case of electron-nucleon scattering. We show explicitly in 
this Letter that the reduced ^*p total cross sections manifest the Froissart In 2 s growth 
with hadronic energy s that also characterizes the jp, n^p and pp and pp cross sections at 
very high energy Our fits to 28 different Q 2 data sets of F^x^Q 2 ) demonstrate that a 
In 2 s growth holds for each of the virtual photon's mass Q 2 . This observation allows us to 
obtain a new constraint on the Bjorken x dependence of the quark and antiquark parton 
distribution functions at small x, viz., they should behave at very small x as (3 \og 2 (xo/x), for 
all Q 2 3> m 2 , where m is the proton mass, and (3 and xq are functions of Q 2 . As x — > and 
Q 2 — ► oo, (3 and x approach constant values /?', x' Q . Only 6 parameters are needed to fit 28 
data sets, allowing us in principle to extrapolate to very large energies and very low x, well 
beyond the experimental range presently accessible, with confidence that we understand the 
x dependence of structure functions and parton distribution functions for "wee" partons. 

Kinematics. In the inclusive process ep — > e'X, the laboratory four- vector momentum 
of the exchanged virtual photon 7* is q = {v,q), with negative q 2 = — Q 2 , and the proton 
laboratory four- vector momentum is p = (m, 0). The invariant s, the square of the center 
of mass (cm.) energy W of the j*p system, is s = W 2 = (q + p) 2 = 2mv — Q 2 + m 2 . The 
Lorentz invariant variables x and y are defined as x = -Q— = ^— and y = where k is 

" 2p-q zmu a p-k ' 

the incoming electron's four-vector momentum. Thus, 
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In terms of x, y and Q 2 , the cross section for ep — * e'X may be expressed as 
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Here, F^x^Q 2 ) is the dimensionless proton structure function, R = f ^ Fl , with being 
the longitudinal structure function, a the fine structure constant, and T = a(2 — 2y + 
y 2 )/(2nQ 2 x). 

The quantity measured experimentally, a^ p (x,y,Q 2 ), is a function of 3 variables, x, y 
and Q 2 . We define the total virtual photon cross section (x, Q 2 ), a function of only two 
variables, x and Q 2 , as 

a\% = a T (x,Q 2 ) + a L (x,Q 2 ) (3) 



7 P 

4ir 2 a fl + 4m 2 x 2 /Q 2 ^ 



F?(x,Q 2 ) (4) 



Q 2 V 1 - x 
»^ X ^), Q*»4 m V. (5) 

In Eq. Q, cr T (x, Q 2 ) and <tl(x, Q 2 ) denote the transverse and longitudinal virtual photon 
cross sections. 

Saturation of the Froissart bound. High energy cross sections for hadron-hadron scattering 
must be bounded by a ~ In 2 s, where s is the square of the cm. energy. This fundamental 
result is derived from unitarity and analyticity by Froissart 2]. Saturation of the Froissart 
bound refers to an energy dependence of the total cross section rising no more rapidly than 
ln 2 s. 

It has been shown that the Froissart bound is saturated at high energies jsj in jp, ir^p 
and pp and pp scattering, and, as we now will show, also in 7*p scattering. 

We choose to work in terms of a dimensionless "reduced" 7*p cross section, a^?* p {yV, Q 2 ) / k, 
where k = 4iT 2 a/Q 2 . Using the parameterization of Block and Cahn 0, we write the reduced 
cross section 1 as 

<*U W, Q 2 )/k = A + f3\n 2 - + cs-°- 5 . (6) 

-So 

The 4 coefficients A, j3, s , and c are functions of Q 2 . We present fits to 29 data sets 
published by the ZEUS [5] collaboration 2 , for Q 2 = 0.11, 0.20, 0.25, 0.65, 2.7, 3.5, 4.5, 6.5, 
8.5, 10, 12, 15, 18, 22, 27, 35, 45, 60, 70, 90, 120, 150, 200, 250, 350, 450, 650, 800, and 



1 The parameterization of Eq. @ is the same as that used by Block and Halzen Q for their successful fits 
of hadronic cross sections for jp, n^p and pp and pp scattering, except for a transformation of variables. 

2 In order to avoid possible normalization differences, we have not included HI data |6( in our analysis, 
although preliminary examination of the combined data sets leads us to identical physics conclusions. 
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1200 GeV 2 . Prior to displaying and discussing our fits, we make three observations about 
the data: 



• An examination of plots of F 2 (x, Q 2 ) versus x for different values of Q 2 shows that 
all data sets are compatible with going through a common point at x ~ 0.09 and 
F 2 ~ 0.41. We call this common intersection the "scaling" point x-p and examine its 
significance below. 

• Inspection of Ff (x, Q 2 ) versus x shows that for x < xp and Q 2 >, 350 GeV 2 , the data 
overlap each other, i.e., within experimental uncertainties, they scale in Bjorken x (see 
footnote EJ). 

• When plotted as oi^/zc versus W, with a logarithmic W axis, the data at all Q 2 
show similar parabolic shapes for x < xp, rising with W, with the shape parameters 
depending slowly on Q 2 , with a negligible inverse-power term proportional to c. 

Introducing x from Eq. ((TJ), defining the parameter xq{Q 2 ) = Q 2 /so, and setting c = 0, we 
find that Eq. © becomes, for x < xp, 

1 — X XqIJI 2 



o-l%(x,Q 2 )/K=A + (3\n 2 



x - 



X 



+ 



Q 2 



(7) 



For Q 2 ^> m 2 and x < xp, we obtain the closed forms 



a"%{x,Q 2 )/K = A + p\n 2 



F 2 p (x : Q 2 



1 — x 



x - 



x 



l-x) [A + f3W 



1 — x 



Xq- 



X 



(8) 
(9) 



where A, (3, and x are functions of Q 2 . 

In order to bolster the above observations quantitatively, we made x 2 fits of Eq. to 29 
sets of the Zeus data, i.e., using the In 2 s parameterization of Eq. © with c = 0, with the 
added requirement that each curve go through its appropriate (Wp(Q 2 ), cr*°! JWp, Q 2 )/k) 



scaling point, where 



Wp(Q 2 ) = W(Q 2 , x P ) =\ —{l~ xp) + m 2 . 

V x p 



(10) 



We fit only those data with x < xp, corresponding to W > Wp, which is our definition of 
high energy data. We treat the 5 highest Q 2 data sets, Q 2 = 350,450,650, 



4 



800 and 1200 GeV 2 , as one set, since their F%(x) data points overlap. From fits to individual 
sets of points with common Q 2 , we find that over the range 3 0.11 < Q 2 < 1200 GeV 2 , and for 
energies W > Wp, each data set with a common Q 2 is fit satisfactorily with the 3 parameters 
A, (3 and sq, with A constrained by 

W 2 

A = a%(Wp,Q 2 )/K-(3\n 2V -^. (11) 

s o 

Our results are illustrated in Fig. ^ where we plot the reduced cross section vs. the cm. 
energy W for the Zeus |5] high energy data, together with our constrained fits. We use 
their published values of F% (x) and for errors, their statistical and systematic errors taken 
in quadrature. All of the data 4 agree well with a In 2 s parameterization over this wide Q 2 
range. Also shown in Fig. ^as large bold points are the 11 different values of W(xp) that 
correspond to the scaling point xp = 0.09, F$(xp) = 0.41. 

The agreement with experiment, over the broad range of Q 2 investigated here, is seen 
perhaps more clearly in the plot of F%(x) vs. x, shown in Fig. El This figure shows the 
same Q 2 sets as in Fig. ^ along with the same 11 constrained In 2 s fits of the form specified 
in Eq. (J7J), with A constrained by Eq. ([lip. Careful examination of Fig. El shows that the 
data are in excellent agreement with the hypothesis that all of the widely different Q 2 fits 
intersect at the point xp = 0.09, F%{xp) = 0.41, i.e., it's is a very good scaling point, being 
satisfied by all of the ZEUS js| sets as well as the BCDMS collaboration points []]]. 

Having a scaling point gives us a universal anchor point — an analyticity constraint 
for the fits to the different sets of values of Q 2 . As shown in Eq. ([lljh only 2 parameters, 
(3 and Xq = Q 2 /sq, are now required for the fit to each Q 2 data set. This reduction in 
the number of parameters constrains their values, making the fit uncertainties considerably 
smaller, e.g., the fractional error in F% due to parameter uncertainties is ~ 5% , for both 
x = 10~ 6 , Q 2 = 3.5 GeV 2 and x = 10" 3 , Q 2 = 1200 GeV 2 . 

Although not shown explicitly here, we find that the coefficients (3 and xq are smoothly 



3 For Q 2 > 1200 GeV 2 , there are essentially no ZEUS data with W > W P (1200 GeV 2 ), so we stop our fits 
at Q 2 = 1200 GeV 2 . 

4 Although not shown in Fig. 2] to avoid loss of clarity, we have also made constrained fits for Q 2 = 2.7, 
4.5, 8.5, 12, 15, 18, 27, 35, 45, 60, 120, 150, 200 and 250 GeV 2 , with similar results, a total of 25 fits 
to independent data sets at different Q 2 . We found a total x 2 — 113.56 for 164 degrees of freedom, 
corresponding to x 2 /d.f. = 0.692. For the 5 combined data sets at the highest Q 2 , we found x 2 = 26.99 
for 27 degrees of freedom, a posteriori justifying our combining them into a single set labeled 1200. 
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FIG. 1: The reduced cross section a^ p (W,Q )/k vs. W, in GeV, with n = 4ir 2 a/Q 2 . The high 
energy Zeus ^| data plotted here have x < xp. The curves are In 2 s fits, a}?* p /K = A+ (5\v? {W 2 / sq) , 
forced to go through the scaling point xp, F^xp). The 11 large symbols, Wp(Q 2 ), correspond to 
the common intersection (scaling) point shown in Fig. |2J The very tiny points, data from the 
BCDMS collaboration pj], are not used in the fits. The data labeled "true Q 2 =" are fit as if they 
had been part of the Q 2 = 1200 GeV 2 set. 

varying functions of Q 2 , both being parameterized adequately as 

\nx = lnx^ + ^p (12) 

h 
(Q 2 ) 

where the constants (3 ' and x' are given by 



P = P'-u^r^ (13) 



P',af = lim (3(Q 2 ),x (Q 2 ). (14) 

The 6 constants (3 ', x' , ai, a 2 , b\ and 62 are parameters that can be obtained in a global fit. In 
a truly global fit, in addition to these 6 constants we need the two scaling point coordinates 
xp, -Ff^p), i.e., a total of 8 parameters, in order to fit all the available data on F^ix, Q 2 ). 
This Letter demonstrates the feasibility of such a fit. 

Attention may be drawn to the fits for Q 2 = 0.11, 0.20, 0.25 and 0.65 GeV 2 in Figs. [Hand El 
These curves do not show different shapes in W — a j k space from any of their counterparts 




FIG. 2: Fits to the proton structure function data, F^{x, Q 2 ) vs. x, for 15 values of Q 2 . The data 
are those of Fig. ^ The curves are the In 2 s fits of Fig. ^ converted to the Fn — x plane. The 



a 



tiny points are BCDMS [7( data. The vertical and horizontal straight lines intersect at the scaling 
point x P = 0.09,F 2 p (xp) = 0.41. The data for Q 2 = 350,450,650,800 and 1200 GeV 2 are all fit 
together with a single curve. 

at larger Q 2 — since they also are parabolic in shape — but they exhibit more structure than 
the other curves when plotted in — x space. This shape difference is explained by their 
very low values of Q 2 . The term 1 + 4m 2 x 2 /Q 2 in Eq. (J1J) becomes substantially greater than 
unity as x — > xp, and in the relation ^-(1—x) = W 2 — m 2 of Eq. ([Tjl. ^-(1— x) ~ 2m(W — m) 
for small Q 2 and x <, Xp; these terms strongly influence the transformation from <t*°p/k — W 
space to .Ff — x space when Q 2 is very small. We find that one closed form fits all of the Q 2 
data for x < xp, i.e., 

1 — x 



F p 2 {x,Q< 



1 + Am 2 x 2 /Q 2 



x 



A + pln 2 



1 — x ( m 
x I 1 + 



X 



X 



Q 2 l-x 



(15) 
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Implications for parton distribution functions. The structure function F 2 (x, Q 2 ) has a simple 
interpretation in the parton model in which the scattering from the proton is due entirely 
to the scattering from its individual constituents. Only the quarks and antiquarks couple 
directly to the electroweak current carried by the virtual photon 7* in DIS. The DIS cross 
section may then be expressed in terms of the probabilities qf/h{x) of finding a quark, 
and qf/h(x) of finding an antiquark, of flavor / and fractional momentum x in hadron h, 
times the cross section for the elastic scattering of that parton. As a consequence of gluon 
radiation in quantum chromodynamics (QCD), the quark, antiquark, and gluon densities 
become functions of Q 2 as well as x, e.g., qf/h{%) —> 1f/h(x,Q 2 ). 

The precise expression for F 2 (x,Q 2 ) in terms of parton distributions depends on the 
choice of the factorization scheme in QCD. The most commonly used schemes are the DIS 
and the modified minimal subtraction MS schemes. In the DIS scheme, all higher-order 
contributions to the structure functions F 2 (x, Q 2 ) are absorbed into the distributions of the 
quarks and antiquarks [9], and we may write, to all orders in QCD, 

F 2 (x, Q 2 ) = £ e 2 x(q™ s h (x, Q 2 ) + $} s h (x, Q 2 )), (16) 
/ 

where e/ is the fractional electric charge of flavor /. The superscripts DIS indicate that 
the distributions are those defined in the DIS scheme. In QCD, gluon radiation removes 
momentum at large x, decreasing the value of F 2 (x, Q 2 ) as Q 2 grows, and builds up the 
quark and antiquark distributions at small x, leading to the qualitative expectation of the 
scaling point that we observe in x at which F 2 (xp, Q 2 ) is independent of Q 2 . 

The Froissart ln 2 (s/s ) energy dependence of the data allows us to conclude, for x < xp 
and Q 2 » m 2 , that 

1 — x' 



xq?$(x,Q') ~ x [A + (3\n 

Moreover, within experimental errors, we find that 

xq^{x,Q 2 > 350 Gev 2 ) « xqffi 



x - 



x 



X 



(17) 



, 1 — x 

X 1 



"0 

X 



(18) 



~ (l-x) x i^A' + P'ln 2 

with Q 2 -independent values A' = 0.40,(3' = 0.050 ± 0.008, and x' = 0.28 ± 0.01. 

Parton distribution functions of quarks, antiquarks, and gluons are generally extracted 
from a global fit |l]] to a wide class of hard scattering data including deep-inelastic scatter- 
ing cross sections. Typically one begins with assumed parameterizations of the x dependence 



of the distributions at a fixed low value Qo and then uses the perturbative evolution equa- 
tions [l(3] of QCD to obtain the x and Q 2 dependences of these distributions for all Q > Qo- 
The non-perturbative form that is assumed in the CTEQ parameterization, for example, 
for these distributions has power behavior at small x, xf(x,Qo) ~ x A . We contrast this 
power behavior with the logarithmic behavior that is shown by the data, suggesting that 
it is productive to redo a global fit program using logarithmic behavior at small x. Over 
a small range of x one cannot distinguish a logarithmic and a power expression with small 
fractional power (e.g., ~ 0.25), but the behavior of these expressions is clearly different when 
extrapolated over a very great range. 

In conclusion, we demonstrated that we can make simultaneous In 2 s fits saturating the 
Froissart bound, to 218 ZEUS datum points, with all fits going through the same scaling 
point, strongly constraining the behavior of F%(x, Q 2 ) at tiny x. These fits are made with 
only 2 parameters per set of data with common Q 2 , since the analyticity constraint ^] of 
Eq. (jH| on the cross section at the scaling point requires only the 2 coefficients, (3{Q 2 ), 
xq(Q 2 ). This additional constraint insures a well determined set of fit parameters, and thus 
an accurate forecast of the very small x behavior. Saturation of the Froissart bound allows 
us to project to very large W and hence to very small x with a high degree of confidence 
in our functional form. Our demonstration of Froissart energy dependence at small x for 
each Q 2 should shed light on various theoretical efforts in examining small x physics in QCD 
and on parton saturation at high density (l^ . Our program for the future includes a global 
constrained ln 2 s fit, as described above, to all available F$ data, in order to obtain the 
necessary parameters and their uncertainties needed to extract accurate parton distribution 
functions at small x. 
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